Astronomy & Astrophysics manuscript no. finalrev 


February 5, 2008 


(DOI: will be inserted by hand later) 





Abundance anomalies in hot liorizontal branch stars of the 
Galactic globular cluster NGC 2808. ^ 

G. Pace^ A. Recio-Blanco^, G. Piotto^ and Y. Momany^ 



o 
o 

(N 



(N 
> 

(N 

\o 

O 

o 

9^ 
6 



X 



* Dipartimento di Astronomia, Universita di Padova, Vicolo deH'Osservatorio 2, 1-35122 Padova, Italy 
2 Dpt. Cassiopee, UMR 6202, Observatoire de la Cote d'Azur, B.P. 4229, F-06304 Nice Cedex 4, France 

the date of receipt and acceptance should be inserted later 

Abstract. We present metallicity measurements of 25 stars in the blue horizontal branch of the Galactic globular cluster 
NGC 2808. Our measurements are based on moderate-resolution spectra taken with the multi-object fiber facility FLAMES- 
UVES, mounted on Kueyen at the Very Large Telescope. We confirm that stars hotter than a threshold temperature have 
super-solar abundance, while the cooler ones respect the nominal metallicity of the cluster, i.e. [Fe/H]^-l.l. The threshold 
temperature is estimated to be about 12 000 K, corresponding to the so called j(-jump, and coincides with the sudden departure 
of the cluster horizontal branch from the models. The metallicity increases with temperature for star hotter than the jump, 
confirming the hypothesis that the process responsible for this abrupt metallic enhancement is the levitation due to the strong 
radiation field in absence of a significative convective envelope. A metallicity dependence of the abundance enhancement is also 
suggested, with more metal poor clusters having a higher increase in metal content. The slope in the temperature vs. abundance 
diagram is higher than the errors involved, and the metal content of the cluster plays possibly a role in determining the amplitude 
of the jump (more metal poor clusters show more enhancement after the jump), although other parameters, such as clusters' 
characteristics and even the atomic species involved, may also someway contribute. 



1. Introduction 

Stars in the Horizontal branch (HB) have already passed the 
red giant branch (RGB) phase and are n ow burning helium in 
their core ( iHovle & SchwarzschildllQ's^ . The wide colour dis- 
tribution of the HB, called the HB morphology, is the result 
of large differences in the envelope mass of stars having the 
same core mass, at the same evolutionary stage. HBs in metal- 
rich clusters are usually populated mostly in the red part of 
the RR-Lyrae instability strip, whereas metal-poor ones tend 
to have a n extension on the blue side, as first recognised by 
ISandage & Wallerstein ( 196^). It soon became also clear that 
clusters metallicit y is not the onl y parameter governing the HB 
morphology ( Sandage & Wilde vi. 1967: .van den Bergh..l967.) . 
As a matter of fact, in several globular clusters (GCs) with in- 
termediate and metal-rich content, a consistent or even dom- 
inant fraction of the HB stars is located blueward of the RR- 
Lyrae instability strip, locus to which we refer as BHB for the 
remainder of the paper 

Various issues that have arisen about the HB of GCs from 
the observational point of view are seriously challenging the 
HB standard models: 



■ some clusters have HB blue tails th at extend very far to the 
blue, up to ~ 30 000 K or more (see lCastellani & Castellani 
1^993, for a theoretical explanation); 

■ the stellar distribution of the HB is characterised by dif- 



ferent gaps 
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jumps, i.e. discrepancies between the observed and the 
modeled distribution of stars in the colour-magnitude di- 
agram ( Grundahl et al. 1999; Momanv et al. 2002); 
■ the rotation-velocity distribution also shows peculiarities: 
some HB stars do not rotate as slow as expected for low- 
mass stars and even a few fast rotators (v sin / ~ 40 km/sec) 
have been found in several clu sters redward of the first lu- 
minosity Jump at ~ 11 500 K dPetersonll 1 983t iBehill 1 999l 
[Recio-Blanco et al. 2004) 

Gradually, some scenarios have been proposed in order to 

elucidate these facts which could be all intertwined (see e.g. 
I — — 11 — 1 

Behr 2003a, and reference therein). Still, the unresolved puz- 
zles outnumber the answered questions. 

The HB luminosity jump was first discovere d in M 1 3, in an 
intermediate-band photometric study made bv 'Grun dahl et alJ 
(.1998) . It was seen as a sudden departure of the observed loca- 
tion of the HB st ars in the (M-v)n vs. V diagram from the model 
track. Soon after dGrundahl et al.ll9'99l) . this feature was recog- 
nised to be ubiquitous, occurring at the same colour, i.e. at the 
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same temperature within the errors, in all the surveyed clusters, 
no matter what their metallicity and morphology is, provided, 
of course, that they extend both blueward and redward of the 
jump. At a temperature higher than ~ 20 000 K, the horizon- 
tal branch position on the colour-magnitude diagram becomes 
again consistent with the model track. This photometric jump 
correspond s with that in th e Tpff v s. log g diagram of M 15 dis- 
covered by iMoehler et al.l ( Il995h . The mechanism suggested 
to be most likely responsible for both features is the radiative 
levitation of heavy metals, that increases the atmospheric abun- 
dance of the star. As a result the opacity due to the metallicity 
rises with respect to the hydrogen opacity. Hence there is an en- 
ergy redistribution in which the Balmer jump is partially filled 
in, but the bolometric emission remains unchanged, and it is 
therefore seen as an increase in luminosit y in the u Stromgren 
and U Johnson magnitudes (that is whv .Bedin et al ] ((2000,) find 
the gap also in the U - B vs. U diagram). The occurre nce of 
levitation was already argued bv iMichaud at all ( Il983l) . They 
predicted that, if the outer envelope of hot HB stars is stable 
enough to al low hehum gravi t ationa l settling, which was hy- 
potesised by lOreenstein et al.l (|l967) in order to explain ob- 
served helium underabundances, then there should be no such 

a mechanism like convection able to prevent the metal enhance- 

1^ , .ly- 



ment caused by the radiation. Caloi (1999) used radiative levi- 
tation of metals to explain the presence of a gap at 8 500 K seen 
in seve ral clusters (but not in NGC 2808'). lHui-Bon-Hoa et al.l 
(12000^ made the only attempt to our knowledge to model the 
atmosphere of the HB stars taking into account radiative lev- 
itation and gravitational settling. They managed to reproduce 
the main observational features, such as the photometric jump, 
since the luminosity in each bandpass depends on the total 
opacity and not on the detailed abundance pattern. However 
these authors warn that the element abundances can only be 
reproduced with more sophisticated models. 

The advent of the new generation telescopes made possi- 
ble to check in a more direct way whether radiative levita- 
tion is at work and to make quantitative analysis. Chemical 
abundanc es on blue HB stars have been mor e extensively re - 
visited by'Behr ("2000a"b'), further extended in'Behr'('2003a'b); 
[Moehler et al. ( 2000) : Fabbian et al. (2005 ) . New detailed mea- 
surements of several elements confirmed that blue HB stars 
cooler than 1 1 000-12000 K in general show no deviation from 
the globular cluster abundances derived from red giants. On 
the contrary, stars hotter than 11000-12000 K, depart from 
the general globular cluster abundances. In particular, iron, ti- 
tanium and chromium enrichments to solar or even super-solar 
values are detected. In addition, some metals, such as phospho- 
rus and manganese, display significantly larger enhancements 
than iron. On the other hand, the abundances of magnesium, 
silicon and calcium, both below and above the critical tem- 
perature of 1 1 000 K, are consistent with very little or no en- 
hancement, even though a large, likely real scatter of silicon 
abundance is observed. Finally, underabundances of helium for 
stars hotter than 1 1 000 K are measured. 

Stellar rotation could be tightly connected to the HB chem- 
ical anomalies. In particular, an abrupt change in the rotational 
velocity distribution along the HB at the level of the jump at 
~ 11500 K (.Recio-Blanco et al...2002. .2004: .Behr..2003a.b,) 



can possibly be associated to the onset of radiative levita- 
tion. Angular momentum transfer prevented by a gradient in 
molecular weight (in a picture a la lSills & Pinsonneau lt 2000) 
or removal of angular momentum due to enhanc ed mass loss 
(Recio-Blanco et al. 2002; Vink & CassisilEooi) could be re- 
sponsible for that discontinuity. In any case, the average ro- 
tation rate is surprisingly high, considered that braking mecha- 
nisms should have spun down the progenitors of these HB stars, 
ever since their entrance on the main sequence, and the overall 
observed rotation distribution is far from being understood. 

In summary, the observational results of the recent years are 
unveiling the d etails of a complicated puzzle of physical ^ aram- 



eters (see e.g. 'Moe hler et al.ll2002t iRecio-Blanco et al 2005h 
and open a number of new questions. Is there a dependence 
of the chemical enhancements with the original metallicity of 
the stellar population? How is stellar rotation influencing the 
chemical abundance anomalies or vice versa? The answer to 
these questions will not only help to our fine-tuning of stellar 
evolution. Hot HB stars play an important role in population 
synthesis of extragalactic non resolved systems; as these old 
hot stars can be easily confused with young hot stellar popula- 
tion. For the same reason, they can also afi'ect the modeling 
of the star formation history in dwarf galaxies of the Local 
Group which tries to reproduce their colour-magnitude dia- 
grams ('Tolstoy"l998V New theoretical models with diffusion 
and radiative levitation on HB stars, supported by appropriate 
measurements, are essential to reconstruct the general picture. 

This work is intended to study the radiative levitation in 
the HB of NGC 2808. With that purpose, a sample of 32 stars 
belonging to the HB of this cluster, in the range of effective 
temperatures from 11000 to 16000 K, has been observed. 
NGC 2808, with its relatively high metaflicity, ([Fe/H]=-1.14 
according to Carretta et al 2004), is one of the most pecu- 
liar examples of second-parameter GCs. It h as a very ex - 
tended blue tail, a marked bimodal distribution (lHarrislll974h . 
being the loci both redward and blueward of th e instabiUty 
strip very well populated (see iBedin etaP bOOO' for a recent 
high-precision photometr y). The RR-Lyrae variables are very 
sparse: Clement & HazerJ( ll989l) foun d only two, while 16 of 
them could only recently be found bv lCrowin et alJ ( l2004h by 
using new image-subtraction techniques that allowed the anal- 
ysis of central crowded field. There are two additional signifi- 
cant gaps in flie B HB star distribution at 18 000 and 25 000 K 
(ISosinetallli997l) . Moreover, enhancement and spread of he- 
hum among cluster stars has b een recently suggested as a way 
to explain its HB morphology JP' Antona et alJl2005l) . Finally, 
NGC 2808 is about 20% younger than metal poor clusters of 
our Galaxy dRosenberg et alJl999tlDe~ Angeli 2005). 

Section 2 of this paper presents the observations and the 
data reduction. In Section 3, the atomic data and the deriva- 
tion of the atmospheric parameters are addressed. The results 
are presented in Section 4 and discussed in Section 5 through 
a comparison with previously analysed clusters. The related 
topic of the rotation velocity w ill be treated in a separate work 
JPaladini et alJin preparationh . with a much wider dataset than 
the one used herin, to compensate with a rich statistic the pro- 
jection of the rotational velocity, and with many more stars 
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cooler than the luminosity jump, where most fast rotators have 
been found in other clusters. 

2. Observations and data reduction 

We have observed, in the ESO observing run 072.D-0742(A), a 
total of 32 stars, selected fr om the ground based wide field pho- 
tometry bv lMomanv et al ](2003, 200'i). Due to their faintness 
(their visual magnitude is between 16 and 17.5 mag) and to the 
fact that they were collected within a more general programme 
aimed at the measurement o f radial velocity of many hundreds 
of stars fPio tto et al J 12004V some of the spectra have a lim- 
ited S/N ratio, ranging from few units to a few tens. Therefore, 
not all the observed stars were suitable for a reliable measure- 
ment of their abundance. Our final sample contains 24 stars, 
18 of which are hotter than 12 000 K. The homogeneity of 
the star distribution in a temperature range that extends to over 
16 000 K and the considerable number of the targets, make our 
sample unique among those so far observed in a single cluster 
for the same purpose. 

Targets lie in the low-crowded outskirts of the GC, to avoid 
contamination from other stars. The spectra were obtained in 
service mode in December 2003 and February 2004. The data 
have been collected with the FLAMES fiber facility connected 
to the UVES spectrograph at the focus of the Unit 2, Kueyen, 
of the Very Large Telescope. The resolution was of R=47 000 
and the exposure times ranged between 2 400 s and 3 000 s. The 
spectra have been reduc ed with the FLAMES-UVES pipeline 
jModigUani et alJl2004l) . and then analysed both with MIDAS 
and IDL routines. 

3. Abundance measurements 

The abundance measurements were performed using the pro- 
gram WIDTH3, developed by R. G. Gratton and ad apted by 
D. Fabbian to temperatures up to Teff ^ 20000 K ("Fabbian' 
12004: Fabbian et al. 2005). The procedure establishes the abun- 
dances of chemical species by reproducing the observed equiv- 
alent widths. Once a starting set of values for the effective tem- 
perature, surface gravity, and model metal abundance is de- 
rived, an appropriate spectrum for each star is ob tained from 
the grid of model atmospheres bv lKumczl ('l994') by interpo- 
lating linearly in the temperature and logarithmically to get 
the Rosseland opacity, electronic and gaseous pressure, and 
density. Continuum opacity was obtained taking into account 
all important continuum opacity sources for stars as hot as 
Teff - 20000 K. CoUisional damping constants were com- 
puted with the Unsold formula. The equation of transfer was 
then integrated through the atmosphere at different wavelengths 
along the line profile and theoretical equivalent widths of the 
lines computed and compared with the observed ones (see 
iFabbian et al..2005t for further details). 

3.1. Oscillator strengths, line selection, and solar 
abundances 

Only lines clean from blends were considered in the analy- 
sis. In order to properly identify the spectral lines, we have 



used a line list suitable for our stars. The list assembled 
and used by iFabbian et al.l ( l2005l) was completed, for wave- 
length values l arger than 5 000 A, thanks to the DREAM 
line data base jBiemont et alJ Il999l) . accessible via VALD 
(lKupkaetalJ|l999r and references therein). For the determi- 
nation of the oscillator strengths (^s), laboratory values were 
considered whenever possible. For Fe i lines they w ere taken 
from papers of the O xford group (ISimmo ns & Blackw ellll982t 
Bard & Kock 1994). The gfs for the Feii li i ies we re taken 
from Heise & Kock Cl99Q> and'Biemont et al.l (1 1999^. The so- 
lar abundances were taken from Grevesse & Sauval (1998) . 

3.2. Equivalent Widths 

The use of a FORTRAN program (courtesy of Patrick Francois) 
to automatically detect lines and measure their equivalent 
widths was attempted, but due once again to the rather low S/N 
ratio of the spectra, a more interactive procedure turned out to 
be more effective. Namely, we have manually detected the pres- 
ence of lines and then fitted those lines that were identified (see 
next Section) with a Gaussian. We have then finally measured 
the equivalent width. The two latter steps were performed with 
IDL. 

Due to the quality of the spectra, the EW uncertainties were 
the most important source of errors in the abundance deter- 
minations. This error was evaluated by comparing the equiv- 
alent width measurements of the Feii line at 4923.93 A in the 
five following stars: 37456, 41077, 41388, 46386 and 8072, 
namely, the stars before the jump for which this line was de- 
tected. These stars have almost equal atmospheric parameters, 
with differences just slightly larger than the internal errors. An 
analysis of the line measured in the five spectra, indicates that 
the standard deviation of the equivalent width measurements is 
12.6 mA. We expect this value to marginally overestimate the 
equivalent width errors, sin ce the s t ars Ar e not identical. The 
value calculated through the ICavrel ( 1 1 9 8 8l) formula: 

A{EW) = 1.6 yJ(wdx)/(S/N) 

with w being the full width at half maximum (FWHM) typical 
of the lines, in this case ~0.12 A, and dx being the pixel size, 
about 0.03 A, we get A(EW) - 10 mA(assuming a S/N -10). 
As expected, this value is lower than the standard deviation of 
the line in the five stars considered, but only by a small amount. 
This result impUes that identification of the correct continuum 
level (that is neglected in the Cayrel formula) is not an issue 
here, as expected since the spectra of blue HB stars have very 
few absorption lines. 

Consequently, we calculated the abundance errors, derived 
from an equivalent width uncertainty of 10.4 mA. 

3.3. Atmospheric parameters 

Model atmosphe res appropriate for each star were extracted 
from the grid of iKuruca ( 1 199 4V T^fi values were determined 
from the c olour of the sta r s in th e Johnson U vs. {U - V) 
diagram by iMomanv et alJ (l2004t see Fig. 1 therein) and in 
the a U vs. (U - B) diagram of the same autho rs. The mod- 
els of Zero Age Horizontal Branch (ZAHB) by ICassisi et alJ 
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Table 1. The photometric data adopted to estimate the effective temperatures of the stars and the relative outcomes. Both U-V 
and B-V colours can be used, obtaining the results shown in the respective columns. 
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Fig. 1. The U, ( U-V) diagram of NGC 2808 from 
iMomany et al.V2004') and the ZAHB model from'C assisi et alJ 
(.1999.) . The target stars Are marked as filled squares. 



i 19991) were used. Hereinaft er, we adopt a reddening value of 
E(B - y) = 0.22, from the iHarrisI ( 1 19961) compilation in its 
revised version of 2003. 



To derive the final effective temperature of the programme 
stars, we calculated the mean of the derived temperatures using 
the (U-V) and the (U-B) colours. Table 1 presents the photo- 
metric parameters for all the observed stars in NGC 2808, with 
the corresponding temperature values derived from the colour- 
magnitude diagram of Momany et al, (2004), labelled Teffu_y 
and Teffij_g. Internal errors in the adopted temperatures were 
derived in the following way: we have derived four values of 
the temperature, two for each of the two photometries, respec- 
tively adding and subtracting the errors on the colours. We have 
then computed the half of the difiference between the largest 
and the smallest values obtained. Clearly, this kind of proce- 
dure relies on present ZAHB models without radiative levita- 
tion unable to reproduce the observed horizontal branches be- 
tween Teff ~ 12000 K and 20000 K. This introduce external 
errors that are difficult to evaluate. In particular, the tempera- 
ture determination relies on the criterion adopted to choose the 
point on the model track that represents our sample star, which, 
if enhanced, is far from being on the track. In the present anal- 
ysis we have taken, for each star, the point on the model track 
that has the same colour as the star, according to the fact that 
temperatures are more colour sensitive in this region. Another 
approach is to choose the closest point of the ZAHB to the star. 
On the other hand, one of the objectives of this work is also 
to help the developing of the models. Nevertheless, in order 
to compare the effective temperatures adopted throughout this 
paper with those adopted in other abundance analysis of blue 
HB stars, we plotted (see Fig. |2} the runs of Teff against the 
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479 
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0.0+1.0-0.0 


53926 
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Table 2. Atmospheric parameters used to carry out the abun- 
dance analysis of our sample stars. For the initial guess on the 
metallicity the cluster's iron abundance, i.e. [Fe/H]=-1.14 dex, 
has been used for the stars before the jump, and the solar value 
for the others. These values have been iteratively changed until 
the achievement of the final convergence within 0.01 dex. 



dereddened Johnson B-V colour for both our p rogram stars i n 
NGC 2808 and the stars in NGC 288 analyzed by Behr (2003a). 
NGC 288 was chosen because it has the same metalli citv as 
NGC 2808. The reddenings here adopted are those from lHairisI 
( l2003h compilation. 



* -ir * 



-0,15 -0.10 -0,05 0,00 

(B-V) „ 



Fig. 2. Comparison between the temperatures derived for 
NGC 2808 and those of Behr for NGC 288, chosen because 
of the similar metallicity. 



As in iFabbian et al ] (l2005 l). gravities are not well con- 
strained by our spectra: in fact equilibrium of ionization is sub- 
ject to possible departures from LTE, and the Balmer lines are 
too broad for reliable determinations of their profile from our 
echelle spectra. We have then determined the gravities of our 
target stars from their cor responding eff ective temperatures, us- 
ing the same method that ,lFabbian et alJ ( 12005,.): a mean re lation 
between log T^ff and gravity (log g), from iBeliJ (l2003al) mea- 
surements of blue HB stars in NGC 288, that were expected 
to be very similar to those in NGC 2808 on the basis of the 
colour-magnitude diagram: 



logg = 2.72- log Teff- 7.10 



(1) 



Errors in the gravities have been derived from the errors in 
the temperatures using Equation^ 

Microturbulent velocities ^ might be derived by eliminating 
trends of the derived abundances with expected line strength for 
some given species. However, given the quality of our spectra, 
generally too few lines could be measured, and with too large 
scatter, to significantly constrain microturbulent velocity. For 
the cooler stars, we relied on Behr's analysis of M 13 stars (as 
no cool HB stars have been already analysed in NGC 288). 
From his analysis, we derived the following relation between 
log Teff and the microturbulence velocity: 



^ = -4.7 • log Teff + 20.9 km/s 



(2) 



Uncertainties in these values for the microturbulent veloci- 
ties can be obtained from the scatter of individual values around 
this mean relation: this is ~ 1 km/s. For the warmer stars, which 
have a very stable atmosphere, still enough to let helium sink 
and heavy metals levitate, we adopted microturbulent veloc- 
ity, in agreement with the values derived by fBehjj (l2003ah in 
NGC 288. 

A check on the accuracy of the chosen atmospheric param- 
eters is given in Figure|3l The plot of excitation potential (x) vs. 
[Fe/H] and Equivalent Width(EW) vs. [Fe/H] is shown for all 
the stars with eight or more Fe ii line detections. The number 
of Fe II lines is never bigger than 13, hence we cannot use this 
method, as already stated, to determine directly atmospheric 
parameters. But the fact than in none of the 5 stars analysed a 
significant slope appears, suggests that ^ and \ogg are properly 
chosen. This is actually a check on the ^=0 assumption for the 
stars after the jump, since all of the five stars in Figure|3]are on 
its blue side. 

Finally, metallicities were obtained varying the metal abun- 
dance [A/H] of the model until it was close enough, namely 
within 0.01 dex, to the derived [Fe/H] value. 

Table 2 presents the atmospheric parameters of the anal- 
ysed stars with their corresponding internal errors. Once Teff, 
log^, ^ and atmospheric [Fe/H] were determined, we cal- 
culated the mean abundance and dispersion for each ele- 
ment. Since the abundances for each element depend upon 
the adopted parameters, we recalculated the abundances while 
varying each of these parameters in turn, in order to test the 
stabihty of our results. 
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Fig. 3. Check of the accuracy of the photospheric parameters 
chosen for the chemical analysis. The lines that on the top- 
panel graphics appear crowded are zoomed in the bottom panel. 

4. Results 

Table |5] contains the abundances of iron, titanium, and 
chromium derived for each of our target stars. The va l ues ar e 
relative to the solar abundances of lGrevesse & Sauv af ("1998). 
The scatter of the measurements for the corresponding lines 
gives the quantity cr^, which is added quadratically to the er- 
ror contributions in Teff, g, ^, and [A/H] described in the pre- 
vious Section, to get the final error for each element in each 
target star. The values for the dominant ionization stage (usu- 
ally singly-ionized metals) were used as a best indicator of the 
actual chemical abundance, since these stages are least suscep- 
tible to non-LTE effects. 

Figure ^ shows the iron abundances as a function of the 
eflBCtive temperature derived for our targets in NGC 2808. It 



is clear that a pattern of radiative levitation is present for blue 
HB stars, with iron abundances increasing gradually as temper- 
ature increases along the HB, up to a value of [Fe/H] ^ -i-l.O 
dex. One more remarkable feature, not so clearly observed in 
previously analysed clusters, is the smooth ascending metal- 
licity trend with temperature after the jump: the levitation be- 
comes more and more eflicient at high temperatures, hence 
for stronger radiation fields. Such a trend could only be dis- 
covered because of the richness of our sample and the uni- 
formity of the temperature distribution of the stars in a very 
wide range, that goes from ~ 1 000 K before the jump up to 
over 16 000 K. These qualities are unprecedented for single- 
cluster samples observed in this kind of studies. The trend can 
be either real or result from an incorrect temperature scale. On 
the one hand, there are to date no available model atmospheres 
that take into account radiative levitation. The discrepancy be- 
tween the predictions of the model and the photometric data in- 
creases with the temperature, as clearly noticeable in Figure^ 
Consequently, abundance overestimations could be stronger for 
stars at larger temperatures, thus producing a spurious trend in 
the abundance vs. Teff diagram. On the other hand, the slope 
that we find in Figure 0] is too large to be easily dismissed 
as an effect of systematic errors and it has, on the contrary, a 
straightforward physical explanation: stronger radiation fields 
cause higher enhancements. We will not be able to say the final 
word on this until models for hot HB stars of stars are refined. 

Although the overall trend, whether real or artefact, is un- 
mistakably apparent, the spread of the data points about it is not 
marginal and there two strong outliers: 14923 at retf=15 606 
and 49334 at Teff ~=13 168, which are respectively over 1 dex 
above and almost 1 dex below stars with similar temperatures. 
The reason could reside in the measurement errors, as they are 
both stars with only one line detected. In addition, 49334 has 
also a large error in the temperature evaluation. We also have 
to remind ourselves that for the hottest stars the non-LTE be- 
haviour of the stellar atmosphere is more and more important, 
hence our measurements, based on LTE modelling, less pre- 
cise. This could lead to spurious trends for the few hottest stars 
of the sample, for which one would be tempted to suggest a 
saturation eff'ect. 

Looking at Figure |4] we notice, in addition, that the onset 
of levitation is quite abrupt for NGC 2808: in less than 500 K 
there is a change of more than 1 dex in Fe. 

On the other hand, as it can be seen from Figure 0] the 
iron abundances for the coolest stars are very close to the clus- 
ter metallicity obtained here from analysis of red giant branch 
stars. We obtain an average value of [Fe/H]=-1.16 + 0.08, 
that is compatible with the abundances obtaine d from anal- 
ysis o f red giants: [Fe/H] =-1.14 on the Carretta & Grattoiil 
(119971) scale (which should be cons i stent w ith the present anal- 
ysis, as the code of 'Fabbian et alJ ll2005l) for our metallicity 
measurements is the ext ension of the software used in the 
ICarretta & GrattonI (Il997| paper); and th e slightly lower value 
of [Fe/H] =-1.23 given by 'Kraft & Ivans' ("2003). 

Due to the limited S/N ratio and the presence of cosmics, 
detections of chromium and titanium Unes with a shape smooth 
enough to reliably measure their equivalent width, could only 
be made for 4 stars. 
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Table 3. Result of the abundance analysis of iron, chromium and titanium and relative errors resulting from the 1 cr uncertainties 
in the assumed parameters. In the last column the resulting 1 cr error is indicated. The quantity in the cr;,„gj column is the standard 
deviation when the number of lines detected is larger than 4. When there are only 4 lines or less, the value reported in the same 
column is half the difference between the largest and lowest measured abundances. 



5. Comparison with other clusters 

The measurement of radiative levitation in clusters with differ- 
ent metallicities is necessary to provide the appropriate data for 
the development of the models. NGC 2808 offers the opportu- 
nity of observing this effect in a cl uster with a very particular 
HB morphology and a younger age ('Pe A ngelil2005l) . 

An important cross check of the ubiquitous nature of the 
phenomenon of the levitation comes also from the comparison 
with NGC 288, a cluster with a metal content very similar to 
that of NGC 2808. 

5.1. Metal enhancements and cluster metallicity 

Figure |5j upper panel, shows the iron abundances derived in 
this paper as a function of the effective temperature, together 
with those of 'Behr ('2003a') for M 15, M 13 and NGC 288, and 
those of M 79 from Jabbian et al.. (.2005i) . Error bars are not 



shown in the Figure for a purpose of clearness. The same data 
are also shown in Tabled] where we have written the abundance 
range that iron, chromium and titanium spread in five clusters 
analysed. The clusters in the Table are sorted according to their 
abundance nominal value, from the most metal poor in the top 
to the most metal rich in the bottom. It can be noticed, look- 
ing at the data about iron, that the most metal poor clusters 
tend to have larger abundance enhancements, suggesting the 
idea that all the clusters tend to reach the same value after the 
jump. The abundance spread within each cluster after the jump 
is very high, much more than for the coolest stars, because of 
the trend with temperature, at least in the case of NGC 2808, 
and possibly because of larger errors. As a consequence we do 
not see a narrower abundance distribution of the stars of all 
clusters after the jump. This trend is not confirmed for titanium 
and chromium, which, however, have a poorer statistic. In any 
case, also for iron, no firm conclusions can be stated at the mo- 
ment. 
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Fig. 4. Apparent discontinuity in the [Fe/H] vs. Teff diagram: 
stars hotter than 12000 K have metalhcities more than 1 dex 
higher than the cooler ones. The increasing metallicity trend of 
the metallicity with temperature after the jump is also evident. 



A M79 
O NGC288 
5*? M15 
+ M15 

□ NGC2808 



A M79 

O NGC288 



+ tvi13 

□ NGC2S08 



A □ e Offl 

□□ □ ^ 



A □ n 

+ + + 

o 

<> □□ □ □ 



A M79 
O NGC288 
^ M15 
+ M13 
□ NGC2S08 



8000 9000 10000 11000 12000 13000 14000 



Fig. 6. Same as Fig. 2 but for the titanium abundances 



as discussed in Section 3, not only differences in photometric 
zero point are there, but also the way of computing tempera- 
tures out of the photometry relies on models that do not match 
the data points. It therefore certainly introduces errors, which 
could differ from study to study. We cannot exclude, however, 
the possibility that the onset of radiative levitation takes place 
at slightly different temperatures in different clusters. Except 
for the hottest non-enhanced star in M 15, NGC 2808 has the 
hottest onset of levitation in any of the clusters in Figure [S] 
Once again, whether this is real and, in this case, what other 
properties are related to the onset temperature of the levita- 
tion, are questions that can only be answered with sufficient 
confidence when models that do match the data points become 
available. 



Fig. 5. Comparison of the iron abundances in blue HB 
NGC 2808 stars with the literature values of M 15, M 13, 
NGC 288 (Behr, 2003) and M 79 .Fabbian et al..(.2005.) . 

Another apparent feature of Figure |5] is the overlapping 
of non-enhanced and enhanced region of the various clusters. 
This is likely due to error in temperature determinations, which, 
as we have seen in Section 3, are quite large for a given sam- 
ple studied homogeneously. They must be considerably larger 
when comparing clusters analysed by different authors, since, 



5.2. The case of NGC 2808 

If we now examine the measurements in NGC 2808, it is clear, 
as pointed out in Section 4, that a pattern of radiative levita- 
tion is also present for blue HB stars. Nevertheless, the iron en- 
hancement for stars between ~ 1 1 500 K and 13 000 K, seems 
smaller than what is found in other clusters. Even NGC 288, 
with a very similar metallicity to that of NGC 2808, but with 
only 4 stars analysed, apparently shows a higher influence of 
radiative levitation, right after the jump, than NGC 2808 stars. 
No measurements of titanium or chromium are available for 
NGC 2808 in that range of temperature. 
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Table 4. Comparison of the abundance ranges before and after the jump. The data about M 79 are from Fabbian et al. (2005), 
those about M 13, M 15 and NGC 288 are from Behr (2003a). For the stars after the jump we also show the range of values spread 
by the abundance jumps, i.e. the enhancement of the star abundance with respect to the nominal value of the cluster (A[M/H]). 
In Behr's sample, NGC 288 has no stars before the jump, so the metallicities in the Table are those from RGB stars, by Gratton 
(1987). In the same way, since we could not detect titanium and chromium lines in stars before the jump, we used measurements 
in RGB stars by Carretta (2005). 



On the other hand, the hotter stars analysed in NGC 2808 
present important signs of radiative levitation, arriving to iron 
enhancements slightly smaller than those of M 15 at lower 
temperatures. The enhancements of titanium and chromium 
for stars with Tefi > 1 3 000 K are also very remarkable and 
compatible or even slightly higher than those observed in other 
clusters. 

In brief, NGC 2808 seems to have a progressive metal en- 
hancement with temperature after Teir ~ 12 000 K, spanning a 
larger range of abundances than other clusters, for which there 
is a more abrupt discontinuity at the level of the jump. Is this 
behaviour due to the higher metallicity of NGC 2808, to the 
complexity of its horizontal branch, or, more simply, to the bet- 
ter coverage (in temperature) of the HB in our analysis? In this 
sense, the comparison with NGC 288, having a similar metal 
content but a different HB morphology and age, can be very 
helpful. Unfortunately, up to date only four stars have been 
analysed in NGC 288 and they are not enough to perform an 
adequate comparison of the abundance anomalies in both clus- 
ters. In the future, it would be very interesting to put in relation 
the couple of clusters NGC 288 - NGC 2808, at a metallicity 
of [Fe/H] ~ - 1 . 1 , with the couple formed by M 79 and M 1 3 at 
[Fe/H] ~ -1.4. 

6. Conclusions 

The objective of this investigation was to shed more light on the 
nature of radiative levitation and diffusion of elements among 



horizontal branch stars and its connection with the physical 
properties of the affected stellar population. Globular cluster 
HB stars come from a coeval population with a very small, if 
any, initial abundance spread. As a consequence they offer a 
unique opportunity of tackling this issue. Studies, including the 
present, made so far on the metallicity of HB stars, lack of ap- 
propriate models describing the atmospheres of the enhanced 
stars. Abundance measurements at this stage, however, point 
to important conclusions which are unlikely to be a mere arte- 
fact of the measurement method: the slope in the temperature 
vs. abundance diagram is significantly higher than the errors 
involved and the metal content of the cluster has probably a 
primary role in determining the amplitude of the jump (more 
metal poor clusters show more enhancement among stars af- 
ter the jump). To our knowledge, model atmosphere able to 
account for the radiative levitation are not yet available, even 
though a theoretical work on the phenomenon has been made 
(see Section [Q. If, on one hand, the analysis here presented 
will be considerably improved when more adequate model at- 
mospheres can be adopted, on the other hand the achieve- 
ment proper models require observational clues regarding the 
amount of metal enhancement and its dependence on physical 
parameters and clusters' properties. This paper is also meant to 
contribute to the collecting of such observational clues. 
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